PROCEEDINGS 

°F SCIENCE 



O 



> 
OO 

(N 

(N 
O 
(N 



X: 



Probing Quark Gluon Plasma by Heavy Flavors 



Santosh K Das* and Jan-e Alam 

Theoretical Physics Division 

Variable Energy Cyclotron Centre, 1/AF, Bidhan Nagar, Kolkata 700064, India 
E-mail: 



santosh@vecc . gov .in, jane@vecc . gov . in 



The drag and diffusion coefficients of charm and bottom quarks propagating through quark gluon 
plasma (QGP) have been evaluated within the framework of perturbative Quantum Chromody- 
namics (pQCD). Both radiative and coUisional processes of dissipation are included in evaluating 
these transport coefficients. The dead cone as well as the LPM effects on radiative energy loss of 
heavy quarks have also been considered. The Fokker Planck equation has been solved to study 
the dissipation of heavy quarks momentum in QGP. The nuclear suppression factor, /?aa and the 
elliptic flow v^^ of the non-photonic electrons resulting from the semi-leptonic decays of hadrons 
containing charm and bottom quarks have been evaluated for RHIC and LHC nuclear collision 
conditions. We find that the observed 7?aa and V2 at RHIC can be reproduced simultaneously 
within the pQCD framework. 
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1. Introduction 

The nuclear collisions at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider 
(LHC) energies are aimed at creating a new state of matter - properties of which are governed by 
thermal quarks and gluons, such a phase of matter is called Quark-Gluon Plasma(QGP). The study 
of the properties of QGP is a field of high contemporary interest and the heavy quarks (HQ) namely, 
charm and bottom quarks play crucial roles in this endeavour. Because HQs (i) are produced in the 
early stage of the collisions and hence witness the entire evolution scenario of the matter, (ii) do 
not dictate the bulk properties of the matter (iii) their thermalization time is larger than light quarks 
and gluons and therefore, the propagation of HQs through a QGP may be treated as the interactions 
between equilibrium and non-equilibrium degrees of freedom. The Fokker-Planck equation (FPE) 
provides an appropriate framework for such studies. In the present work we evaluate the nuclear 
suppression factor and elliptic flow of heavy flavours by solving the FPE for the motion of charm 
and bottom quarks in QGP and contrast the theoretical results with the data to extract the properties 
of QGP 

2. Formalism 

The evolution of HQs momentum distribution function while propagating through the QGP 
can be described by the FPE, which reads 



dt dpi 



Mp)f+^[Bij{p)f] 



where the kernels A,- and S,,- are given by 



B 



(2.1) 



Ai = J d^k(o{p,k)ki 

d^k(o{p,k)kikj. 



(2.2) 



The function (o{p,k) is given by 



{271) 



(2.3) 



where gj is the statistical degeneracy, fj is the phase space distribution of the particle j, Vij is the 
relative velocity between the two collision partners and a denotes the cross section. For | p | — 0, 
Ai — jpi and Bjj — D5,y where 7 and D stand for drag and diffusion co-efficients respectively. 

As mentioned before we include both the collisional and radiative processes for HQs dissi- 
pation in QGP. For the collisional processes [^, ^] : gQ — )• gQ and qQ — )• qQ and qQ — )• qQ are 
considered in evaluating the drag (Ycoii) and diffusion (Dcoii) coefficients. 

For the radiative loss the drag coefficient is defined through the relation. 



-dE 
dx 



YradP 



(2.4) 
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and use the Einstein relation Drad = TMy^ad to obtain the diffusion coefficient, where M is the 
HQ mass, T is the temperature of the QGP. In the radiative process the dead cone and Landau- 
Pomerenchuk-Migdal (LPM) effects are included [^. Although the collisional and radiative pro- 
cesses are not independent from each other, however, in the absence of any rigorous method, we 
add them up to obtain the effective drag coefficients, Jeff = Jrad + Jcoiu and similarly the effective 
diffusion coefficient, Deff = Drad +Dcon- This is a good approximation for the present work be- 
cause the radiative loss is large compared to the collisional loss. In evaluating the drag co-efficient 
we have used temperature dependent strong coupling, a, [Q]. The Debye mass, ~ g{T)T is also 
a temperature dependent quantity used here to shield the infrared divergences arising due to the 
exchange of massless gluons. The FPE has been solved with these effective transport coefficients 
and the initial momentum distributions of HQs are taken from the proton-proton (pp) collisions [Q] 
at RHIC and LHC energies. 

The solution of the FPE aie convoluted with the fragmentation functions [Q] of the HQs to 
obtain the pj distribution of the D and B mesons. The transverse momentum spectra of the elec- 
trons originating from the decays: D — Xev and B — )• Xev have been obtained by using standard 
procedures [^. 

(2-1-1) dimensional relativistic hydrodynamics with boost invariance [^] along the longitu- 
dinal direction has been used for the space time description of the of the flowing QGP. We have 
taken the Equation of State (EoS) as P = c^e, where P is the pressure, £ is the energy density 
and Cs is the velocity of sound in QGP. Variation of Raa with i.e. with EoS will be presented 
here. We have assumed the value of the initial temperature, 7] = 400 MeV and thermalization time 
of the QGP, T, = 0.2 fm/c for Au-i-Au collisions at RHIC. The corresponding values for 7] and T, 
for LHC are taken as 700 MeV and 0.08 fm/c respectively. It is expected that the central rapidity 
region of the system formed in nuclear collisions at RHIC and LHC energies is almost net baryon 
density free. Therefore, the equation governing the conservation of net baryon number need not be 
considered here. 

The total amount of energy dissipated by a HQ in the QGP depends on the path length it 
traverses. Each parton traverse different path length which depends on the geometry of the system 
and on the point where its is created. The probability that a parton is produced at a point (r, 0') in 
a QGP of ellipsoidal shape depends on the number of binary collisions at that point which can be 
taken as: 



where R is the nuclear radius and e is the (spatial) eccentricity of the ellipse. A parton created at 
(r, 0') in the transverse plane propagate a distance L = \/R^ — r^sin^<p' — rcos^' in the medium. We 
use the following equation for the geometric average of the integral which appear in the solution of 




(2.5) 



and 



1 




(2.6) 



the FPE [||] involving drag coefficient: 



r = 



J rdrd^'P{r, f) 



dxy{x) 



(2.7) 
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where v is the velocity of the propagating partons. Similar averaging has been performed for the 
diffusion co-efficient. 



3. Results 

In this section we contrast the theoretical results with the experimental data available for 
Au+Au collisions at RHIC energies and also present theoretical results for LHC energy. 

3.1 The Nuclear Suppression 

The transverse momentum distribution of electrons from the heavy flavours produced in p-p 
collisions can be estimated from the charm and beauty quarks distributions which provides initial 
condition to the FPE. The solution of FPE for the HQ transverse momentum distribution contains 
the effects of drag (quenching) on the HQ whereas the initial distributions of HQ does not contain 
any such effects, therefore the ratio of these two quantities, the nuclear suppression factor, Raa 
acts as a marker for the momentum degradation in the medium, which is observed experimentally 
through, Raa{pt) defined as: 

^/ye Au+Au 

RM = /^,„. (3.1) 

The experimental data from both the PHENIX and STAR collaborations [[l^, |ll|] shows substantial 
suppression (Raa < 1) for pr GeV indicating significant dissipation of HQ energy in the QGP. 
The theoretical results describe the data reasonably (Fig. H, left panel). In the right panel we display 
the results for LHC for two values of c^. We have taken an equation of state with velocity of sound 
lower than the value corresponding to the Stefan-Boltzmann limit. Lower value of c., makes the 
expansion of the plasma slower enabling the propagating heavy quarks to spend more time to 
interact in the medium and hence lose more energy before exiting from the plasma. 

3.2 Elliptic Flow 

The elliptic flow, v^^ is defined as: 

V2 {pt) = {cos{2<^)) = ■ ' "^^^ (3.2) 

J dydprdij) \y=0 

In the left panel of Fig. ^ we contrast the experimental data obtained by the PHENIX [11] 
collaborations for Au + Au minimum bias collisions at y^SNW = 200 GeV with theoretical results 
obtain in the present work. It is observed that the v^^ first increases and reaches a maximum of 
about 8% then saturates for pr > 2 GeV. We also find that the data can be reproduced well by 
including both radiative and coUisional loss with = 1/ \/4. In the right panel we display the 
variation of v^^ with pr for LHC energies with Cs = l/V^-. 
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Figure 1: Left panel: variation of Raa with pr obtained from the present work compared with the experi- 
mental data measured by STAR and PHENIX collaboration at ^5n« ~ 200 GeV. The experimental data of 
STAR and PHENIX collaborations are taken from [101 and [VI \ respectively. Right panel: variation of Raa 
with pT for LHC. 




Figure 2: Left panel: variation of with pr at the highest RHIC energy. Experimental data is taken 
from 1 11 1. The value of the "effective" impact parameter, b = 10.2 fm/c for minimum bias Au+Au collision 
at y^Sf^N — 200 GeV. Right panel: variation of v^^ with pr for LHC. 



4. Summary and discussions 

In summary, we have simultaneously reproduced the measured nuclear suppression and elliptic 
flow of heavy flavours at RHIC energies within the framework of pQCD interactions of the non- 
equilibrated HQ with flowing QGP. We observe that the inclusion of the coUisional and radiative 
processes of dissipation and use of non-ideal EoS are the key factors for the successful reproduction 
of the data. Our analysis admits the formation of QGP in Au-i-Au collisions at centre of mass energy, 
^/s^ = 200 GeV with initial temperature ~ 400 MeV. 

Some comments on the R^a vis-a-vis v^^ are in order here. The Raa contains the ratio of pr 
distribution of the electrons resulting from Au-i-Au to p-i-p collisions, while the numerator contains 
the interaction of the HQs with the flowing QGP, such interactions are absent in the denomina- 
tor. Whereas for v^^ both the numerator and the denominator contain the interactions with the 
medium, resulting in some sort of cancellation. Therefore, K factor which is sometimes used for 
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the reproduction of the data may not be same for Raa and v^^ . 

Several theoretical attempts have been made to explain Raa and V2, where the role of hadronic 
matter has been ignored. However, to make the characterization of QGP reliable the role of the 
hadronic phase should be taken into consideration and its contribution must be subtracted out from 
the observables. Recently the drag and diffusion coefficients of hot hadronic medium consist- 
ing of pions, kaons and eta using open charm and beauty mesons as a probe have been evalu- 



ated [|1J, |13[ [14|, [15|, [1^]. It is observed that the magnitude of both the transport coefficients are 
significant, indicating substantial interactions of the heavy mesons with hot hadrons, which may 
have significant consequences of Raa and . 
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